Josephson junction between anisotropic superconductors 
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The sin-Gordon equation for Josephson junctions with arbitrary misaligned anisotropic banks 
is derived. As an application, the problem of Josephson vortices at twin planes of a YBCO-like 
material is considered. It is shown that for an arbitrary orientation of these vortices relative to the 
crystal axes of the banks, the junctions should experience a mechanical torque which is evaluated. 
This torque and its angular dependence may, in principle, be measured in small fields, since the flux 
penetration into twinned crystals begins with nucleation of Josephson vortices at twin planes. 



I. INTRODUCTION 



The high-T c superconductors are often used as poly- 
crystals made of anisotropic grains touching each other 
with various degrees of crystallographic misalignment. 
As a rule, the contacts have Josephson properties. Hence, 
physical characteristics of Josephson junctions with mis- 
aligned anisotropic banks are of considerable interest. 
To our knowledge, the only publications on this subject 
were by one of us and were concerned with a relatively 
simple situation of a junctionJpetween perfectly aligned 
anisotropic superconductors .BH 

In this paper we consider two misaligned grains of the 
same uniaxial material forming a Josephson contact. The 
material is characterized by the dimensionless "mass ten- 
sor" with eigenvalues m a = mi, < m c which are normal- 
ized so that m 2 a m c = 1. The mass ratios are defined as 
Trii/mk — Af/A|, so that the actual penetration depths 
can be written as A-j = \^fm~i where A = (A^Ac) 1 / 3 . 

The paper is organized as follows. We first consider 
aligned banks not only to establish notation and to 
demonstrate the approach; we have in mind a material of 
YBCO family which usually contains sets of nearly par- 
allel twinning planes having Josephson properties. We 
show that the gauge invariant phase difference still sat- 
isfies the sine-Gordon equation, however, the squared 
Josephson length acquires tensor properties which are 
discussed in detail. Next, the case of a general misalign- 
ment is treated and the tensor for the squared Josephson 
length is expressed in terms of mass tensors of the banks 
and their misalignment. 

In twinned materials of the YBCO kind, in increas- 
ing from zero magnetic field, the first vortices entering 
material are usually situated on the twin planes and 
have Josephson properties. We consider the problem 
of a Josephson vortex between anisotropic grains and 
evaluate the torque which should act on the system of 
such vortices and via them on the sample. Although this 



torque is small, given recent developments in improving 
sensitivity of the torque magnetometryjj it is a measur- 
able quantity, from which the Josephson characteristics 
of grain boundaries can, in principle, be extracted. 



II. MAIN EQUATIONS 
A. Aligned banks 

Let us first consider the simplest possible situation: a 
Josephson junction with superconducting banks made of 
the same anisotropic material; rriik are the same on both 
sides and the banks are perfectly aligned. Moreover, let 
the axis y perpendicular to the junction plane xz be the 
principal direction, say b. Let us assume further that a 
Josephson vortex is directed at an angle 9 relative to the 
crystal axis c. We choose the z axis along the vortex. 

The coordinates x, z are related to a,c: x = a cos 9 + 
csin# and z = —a sin 9 + ccosO. Then the nonzero com- 
ponents of the mass tensor in the frame xyz are: 



m„ cos 



m„ sin 



m r sin 



m r cos 



(m c — m a ) cos 6* sin 6* , m yy = m. 



Note a useful property 



m xx m z 



(1) 



(2) 



m 2 a m c . 



which follows from det(mjfc) 

We proceed by calculating the London energy Fl as- 
sociated with a Josephson vortex along z in terms of the 
gauge invariant phase <p. We then minimize the total 
energy Fl + Fj where 



<l>ojc 

2?rc 



dx (1 — cosip) , 



(3) 



with respect to tp to obtain the tunneling current j y in 
terms of if along with the equation governing the phase; 



here cj>o is the flux quantum and j c is the Josephson crit- 
ical current density. 

For a straight vortex along z, all gauge invariant quan- 
tities are z independent: the field components h x ,z are 
functions of x and y, h y — 0, while ip depends only on x. 

We begin with the general expression for the London 
energy 



r 



i 



F -ik^ 



4ttA 2 



mikjiCmlkh)dV 



and integrate the kinetic part by parts to obtain 
„ A 2 



2c 



mikeum f jih m dSi , 



(4) 



(5) 



where dS = ±dxy is an area element on two junction 
sides. This gives 



where the length Xj is the characteristic distance at 
which H changes in the junction plane. In this case 
one can consider the field penetration into the banks 
as the standard Meissner problem of a uniform "ap- 
plied" field H = {H a , H c } penetrating in two half-spaces 
y > and y < 0. Since h a = _ff a exp(— |y|/A c ) and 
h c = H c exp(— \y\/\ a ), we obtain 

h x = H a e-^ IK oasd + H c e~^ IXa sine , 

h z = -H a e- M/X <-- sm9 + H c e- M/K cos9 , (13) 

where A Q , C = \jm a ^ c . In terms of H x%z we have 

h x = H x ( e -\ y \ /K cos 2 9 + e - M/K sin 2 i 

+ H z (e- lyWXa - e ~ M/x A sinfl cos , (14) 



A 2 f°° 
F L = -z- / dx {H x [m zk jk] - H z [m xk j k ]} 

2 C J-oc 



(6) 



h z = H x (e~ lvl/Xa - e~ M/x A sinflcosfl 
+ H z (e-\v\/>* sin 2 9 + e~^ ,Xa cos 2 6 



(15) 



per unit length in the z direction. Hereafter H XjZ (x) are 
the field components at the junction plane y — 0; y is di- 
rected from the side 1 to the side 2 of the junction. The 
symbol [...] denotes the difference of a quantity on two 
sides, e.g., [m xy ] 



The field is continuous at the junction plane y = 
whereas the tangential currents j x<z are not. One finds: 



,(2) 



,(D 



We now utilize the London equations 

4ttA 2 



Ai+ fo dx 

2tt dxi 



-rriikjk 



(7) 



4tt 
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sin 29 



where x is the phase of the order parameter. Also, we 
use the gauge A y — and the continuity of the tangential 
component of the vector potential to obtain 



- 2H Z (— sin 2 9 + — cos 2 O) , 

V A c A„ 

[j 2 ]=2^f^cos 2 

( 



' \A r A a 

// ( — - — ] sin 29 



1 1 
^ A c A a 



(16) 



(17) 



m xx [j x ] +m xz [j z 



ctfio dep 
'8ir 2 X 2 ~dx~ 



in 



:\jx] +m zz [j z ] = 0: 



(8) 
(9) 



Note that the coefficients of H xz in these expressions 
do not reduce to components my, of Eq. (pi) because 
they contain 1/ Jm a ^ instead of m a ^ c . It is convenient 
to introduce another tensor \i a p (a, /3 = x, z) with eigen- 



here we introduced the phase difference <p = \x\, utilized valuef > V y/m^ (Jo obtain M^eplace m ac in Eq. ^ 



the alignment, and the condition [j y ] — 0. Note that 
for the uniform along z solutions, the discontinuities of 
the tangential currents j x and j z at the junction are pro- 
portional to each other. Solving the system (pi), (0) we 
obtain: 



y fj replace 
with l/y/Tn a _ c ). Then, Eqs. (jig), (|17| ) take a compact 
form: 



[J., 



■ [:iz 



c(!)om zz 
T 2 X 2 m a m c 



V(s) 



(10) 



2n 

— \jx\X = H x /i xz — H z /i xx , 
c 

2ix 

— [j z )X = H x fi zz - H z fj, xz . 
c 



(18) 
(19) 



where the property (0) has been used. 

With the help of Eqs. (|8|) and (g) the London energy 
takes the form: 



Ft = 



16tt 2 



dx H z (x) ip'(x) 



(11) 



To make further progress one has to turn to the field and 
current distributions in the banks. Assume that 



Since the ratio [7'x] /[?'*] ^ s nxe d by Eq. (|9[), we obtain a 
fixed ratio H x /H z : 

H x (n xz m xz + \i zz m zz ) = H z (fi xx m xz + fi xz m zz ) (20) 

Clearly, the transverse field H x vanishes if the vortex axis 
coincides with either a or c. 

We now go back to Eq. (p0|) and obtain using Jig) and 



A< A. 



(12) 



W l( \ TJ 

— -~<p (x) = H z , 
AnX 



(21) 



where 



A 



A 



A 



fj, c sin 9 + [i a cos 2 9 li x 



Thus, the London energy 
F L = 



64tt 3 A 



2 /-oo 

[ I dx{y'f. 



(22) 



(23) 



We complement this with the Josephson energy (m and 
minimize the sum with respect to the phase to obtain the 
result of Ref. § 



X 2 jH xx ip" = sin ip, 
where the length Xj is defined as 



(24) 



A, = 



C(j>Q 



16n 2 Xj c 



(25) 



It follows from Eq. (|24|) that the Josephson length is 



where 



/'•. 



Aj(9) = Xj^%j9) , 



:x = k- 2/3 sin 2 9 + k 1/3 cos 2 9, 



(26) 



(27) 



and k = A c /A a is the anisotropy parameter. Note that 
Xj has now meaning of an average Josephson length. We 
see therefore that even in the case of aligned banks con- 
sidered here, the Josephson length Aj(9) depends on the 
orientation of the Josephson vortex within the junction 
plane ac. In particular, we have 



Aj(tt/2) 



Vie. 



(28) 



Using Eqs. (ra) and ( |23j ) and the classical vortex solu- 
tion of Eq. (|2J|), tp(x) — 4tan _1 [exp(a;/A,/)], we calculate 
the line energy of the Josephson vortex 



1 



4 7 r 3 /fc 1 /3 XX.J 



-^-Vsm 2 9 + kcos 2 9. 



(29) 



B. Misaligned banks 

Consider now a simple misalignment: the banks are 
made of identical materials with the axis b on both sides 
along y as before, but the a, c axes are rotated round b 
through different angles on two junction sides. We again 
look for vortex solutions uniform along the vortex direc- 
tion z. The mass tensors on two sides are given by Eq. 
(Ef), however, the angles 9\ ^ 6*2- Equation (0) still holds, 
but in Eqs. (18) and (K5I) the components of m^ cannot be 
taken out of the square brackets denoting the differences 
between quantities enclosed on the two sides: 



[1TlxiJi\ 

{m Z iji] = : 



apQ dip 
~8Tr 2 X 2 ~dx~ 



(30) 

(31) 



note that the summation index i = x,y, z can be replaced 
with the two dimensional (2D) a — x, z because m ay = 
on both sides. We now see that Eq. (pTl) holds, too. 



The field distribution in the bank y > is obtained 
replacing \y\ and 9 in Eqs. (|14|), (p~5[) with y and 6*2, 
whereas for y < 0, \y\ — > — y and 9 — ► 9\. Then, one 
obtains for the currents at the junction side y = +0: 



4ttAj+ = c(H x ntz ~ H zf4x 
iirXjf = c(H x ^+ z - H z jj+ z 



(32) 
(33) 



These relations can be written in a compact form with 
the help of the two-dimensional unit antisymmetric ten- 
sor e a(3 : 

4irXj+/c=-n+ e fSj H 1 . (34) 

For the other side of the junction at y = —0 we have 

4ir\j-/c = fi~ /) efaH 7 . (35) 

To proceed, we rewrite the system (J30|), ( pi]) as 

c4> a 



lapjp] 



T ip'(x)S x 



8tt 2 X 2 
substitute here Eq. (|34|), (|35|), and obtain 

Vaijtiji = -TT^P (X)0 xa . 



Here 



Pap — {m a7 ^ y s}es/3 



(36) 



(37) 



(38) 



and {...} denotes the sum of a quantity enclosed taken 
on two junction sides. 

The system of two equations ([57]) can be solved for 
H x . z : 



rr _ 00 Pzx ,, , 

z ~ 2nX det( Pa0 f (X) ' 
This coincides with Eq. (eTJ) in which 

det(p a/3 ) 



A=-A- 



^Pzx 



(39) 



(40) 



To have the following equations in a more symmetric 
form, we introduce a tensor 



Qfiis ^fiaPau 



1a 1 ^a/3M/37 J ^"yu 



(41) 



instead of the pseudotensor p; a compelling reason for 
this is given below. Then 



A = A 



det(q a 



P) 



^Qxx 



(42) 



and proceeding as above we obtain the sine-Gordon equa- 
tion for the phase: 



. , d 2 ip 
A^^-rrr — 81X1 (f , 



i 2 

xx dx 2 



a 2 _ \2 ^Qxx 



XX " J 



det(q af j) 



(43) 



For aligned banks (6\ = 6 2 ) this coincides with Eq. 
(B2n. To see this, note that in this case m and jl are the 
same on both sides and {m a ^^,p~ ( \ — 2?n Q ^/i / g 7 . Then, 
p zx = —2m zl p lz = — 2(m a /x sin 9 + ra c [i c cos 2 9) and 
det(p a p) = 4m a /i a m c /i c (choose the principal axes to cal- 
culate this invariant). 



C. General equation for the phase 

We derive now an equation for the phase difference 
ip(x, z) which holds at the junction plane xz with no ref- 
erence to a particular vortex solution. Since we no longer 
have uniformity in the z direction, Eq. (O) takes the form 

A 2 f 
F L = — dxdz (H x [m z kjk] - H z [m xk j k ]) . (44) 



Using the London relations (0) we obtain instead of (36) 

afio dp 



[m a i3jpl 



8tt 2 A 2 dx a 



(45) 



The current distributions are still given by Eqs. (34) 
(pq); the system ( p7\ ) now reads: 

0o dip 



Vxx^x ~r Pxz-Hz 

PzxH x + p zz H z = 



2w\ dx ' 

0o dp 
2ttA dz 



(46) 

(47) 



with the tensor p a p defined in Eq. (|38[ ). Solve this for 
H xz and substitute to the energy 



F, = - 



dxdz 



( _ ^^f_ 

[PzZ Pxx) o ~T 

OX oz 



Pzx \dxJ 



(48) 



327r 3 Adet(p a/ a) 

Since the scalar of energy cannot depend on the choice of 
coordinates, the coefficients by the derivatives of p must 
form a tensor. This is the tensor q a p introduced above 
in Eq. @: p xx = q zx , p xz = q zz , p zx = -q xx , and 
Pzz = —qxz- We then have 

dp dp 



Fj = 



dx dz q a f : 



(49) 



327r 3 Adet(g Q/ 3) J W>J dx a dxp 

Minimizing with respect to p the sum of this energy 
with 



Fj = 

one obtains 

d 2 p 



t>o3c 
2?rc 



dx dz (1 — cos ip) , 



A 



Q/3 dx a dxp 



snip . 



A 2 - \ 2 



2q a /3 
det(q a/3 ) 



(50) 



(51) 



III. JOSEPHSON VORTICES IN FLAT SAMPLES 



A. Infinite slab 



Consider an infinite thick slab with z axis perpendicu- 
lar to the slab surface and with a set of parallel Josephson 
contacts at planes y = 0, ±d, ±2d, ... Clearly, the system 
for which this model is relevant is a platelet of YBCO 
with a set of parallel twins. Let the magnetic field 7i be 
applied in the plane xz of the junctions at an angle a to 
the z axis. The fields of our interest are small enough so 
that only Josephson vortices are nucleated at the junc- 
tions while Abrikosov vortices are absent. For a general 
tilt a of the applied field, the vortices are tilted as well 
at an angle (3 still to be determined. We denote as N 
the line density of Josephson vortices at each junction so 
that the average magnetic induction is 



*-$■ 



(52) 



The macroscopic boundary conditions for the slab ge- 
ometry are 



B z —Hz , H x = H x . 



We utilize the first of these to obtain: 



N 



(53) 



H z = Bcos/3= — 4> cos(3, (54) 



whereas 



B x = B sin j3 = 7i z tan j3 . 



(55) 



Since we are interested in the limit of vanishing den- 
sity N of vortices, we can disregard their interaction and 
write the Helmholtz free energy as 



F = ±M 



(56) 



where the line energy e of a Josephson vortex is given in 
Eq. (pi). With the help of Eq. (pi) we write: 



(57) 



where 



F = —H z H Jk + t&n 2 {3. 

4-7T 



Hn = 



TT 2 k 1 / 3 XX, 



(58) 



We now turn to the thermodynamic potential which is 
minimum in equilibrium at a given TL. Given the bound- 
ary conditions (p3T), this potential is 



G — F — B^H'r . 

An 



Since 5F = H ■ SB /An, we have 



5G = —(H Z SB Z - BJH X ) , 

47T 



(59) 



(60) 



i.e., the potential G is minimum indeed at fixed B Z ,H X 
or Ti z , Ti x . It is worth noting that in fact G is the Gibbs 
energy with a subtracted term B Z H Z /An = H 2 /An which 
is constant for the slab geometry. 

Using Eqs. (|59|), (p7|) , and ( |53| ) we obtain 



G(H, 0) = '-P- {H \/k + tan 2 (i-H x tan/3) . (61) 

An V 

Minimization with respect to tan/3 yields the direction 
of vortices: 



tan 2 /3 = 



kUl 



H 2 - HI ■ 



(62) 



We observe that when the direction of the applied field 
changes from a = to a = tt/2, the vortices (B) rotate 
from /3 = to /3 m ax such that 



tan 2 /?„ 



fcW 2 



ff 2 -^ 2 



(63) 



It is of interest to see how the equilibrium vortex den- 
sity changes during rotation of the constant in value ex- 
ternal field H: 

B*=Hlt^a + Hl = Hl H2 ° t? 1 ^; 1] ■ (64) 

n n z 

We see that when a increases, B increases, goes over a 
maximum, and tends to zero for a — ► n/2 {f3 reaches 

PmaxJ ■ 

The Gibbs potential G(H X ,'H Z ) is now readily found: 



H cos a \ / H% — TL 2 sin 2 a . 



Vk, 

An 



In particular, the torque r y = —dG/da reads: 
Vk „ , . H$ + H 2 cos 2a 



t„ = TL sin a 

tf Air 



H 2 - H 2 sin 2 a 



(65) 



(66) 



Note an unusual feature: when a — ► 7r/2, the torque goes 
to a constant value: 



rOr/2) = ^ v^T^. 



(67) 



For a > tt/2 (or 7Y Z < 0), one has to change sign in 
Eq. (p7|) or to replace H z with \H Z \. One can follow 
the same derivation to see that \H Z \ appears in Eq. ( |65| ) 
for the Gibbs potential. In other words, the potential 
G(a) ex \ir/2 — a\ in the small vicinity of a — tt/2. A 
similar situation has been discussed in Ref. |^. 

The nonanalyticity of G and r at a — tt/2 is an artifact 
of the infinite slab geometry. To treat better the vicinity 
of a = 7r/2 we turn to a more realistic sample shape. 



B. Oblate ellipsoid 

The thermodynamic potential which is minimum in 
equilibrium for this case is given in Ref. 0: 



FT Ft \ - 

F = F -M-H 

An 2 



(68) 



All macroscopic fields inside the ellipsoid are uniform 
and related to the applied field by 



Hi + n ik {B k -H k ) 



Hi 



(69) 



where n^ is the demagnetization tensor. Let us consider 
the sample as an oblate ellipsoid with z = c as the axis 
of rotation. Then n z = 1 — 2n x and 



(1 - n)H x + nB x = H x , 
2nH z + (1 - 2n)B z = H z , 



(70) 
(71) 



where the subscript x at n x is omitted for brevity. The 
free energy in this case reads: 



00 

Ho 



BHo 

An 



sin 2 (3 + k cos 2 j3 ■■ 



= ^VBTTkBl. 

The field Hi = An dF/dBi is now readily found: 

Br. 

H x = H 

H z = Ho 



VB 2 x + kBf 

kB z 
VBj+kB 2 ' 



(72) 

(73) 
(74) 



In principle, Eqs. ©, @, @, @ determine both H 
and B in terms of the applied field. This determination, 
however, is cumbersome in the general case. 

To establish at what field H for a given orientation 
a, the vortices start to penetrate the ellipsoid, we first 
exclude Bi from the system (73), (|74]): 



H 2 C + H 2 /k = H 2 . (75) 

At the Meissner boundary, Eqs. (^0h and (FrO) yield 



(1 - n)H x 
Therefore, 

H r 



h x 



2nH z = U z 



HI 



(1 - n) 2 An 2 k 



H 2 



sin a 



Arflk 



H li 



= m 



(76) 



(77) 



(78) 



In the following we restrict ourselves to the case of a 
flat sample, n< 1, and to a narrow angular domain near 
a = 7r/2. It is in this domain, the torque evaluated for 
an infinite slab (n x = 0) has a discontinuous jump. The 
simplification here comes about because H z <§C Ti x along 
with B z <C B x and H z <C H x . Using this and utilizing 
the smallness of n, we obtain from Eqs. (|70|), (fri|), (73), 



B x S= 


H x - 


-H Q (1- 


-n) 


B z 




n 


J 


H x ? 


a H , 






H, 



n z 



l-2n' 



H z n 



i?ofc 



Wa; — H 



(79) 
(80) 



Note that according to Eq. (|78|), for a close to 7r/2, the 
field of first penetration is Ho[l — n) and therefore B x 
given above is positive. 

It is now a matter of a simple algebra to express the 
potential F of Eq. (1681) in terms of the applied field: 



F 



H X (H ~ TL X ) 
87m 



(81) 



In terms of the angle 7 = a — 7r/2 between the applied 
field and the ab plane, we have 



F 



7i cos 7(7^0 ~ H cos 7) 
87m 



This yields the torque for 7 <^ 1 
dF n 



(h, 



8nn 



{2H-H Q ) 1 . 



(82) 



(83) 



Thus, as expected, the torque is continuous at 7 = 
(a = 7r/2) where the sample is in the stable equilibrium. 
The torque is fast increasing in magnitude when |7| in- 
creases. 

The crossover between the regime described by Eqs. 
(pq) and (J67J) takes place around 7,„ which can be roughly 



estimated by equating the torque given in Eq. (p7^ and 
that of Eq. (||): 



Jr, 



2nVk. 



The maximum torque then is 



4it 



Vk 



47r 3 AAj 



H. 



(84) 



(85) 



For A ~ 10 _5 cm and Xj ~ 10~ 4 cm, r m ~ 10 2 erg/cm 3 
in fields on the order 100 G. Even for a tiny crystal with 
dimensions (0.1 x 0.1 x 0.01) mm 3 = 10~ 7 cm 3 , we es- 
timate T m as 10 -5 erg. The sensitivity of piezoresistivc 
torque magnetometers is in the range of 10 -7 erg,B so 
that the torque we have calculated here can, in principle, 
be measured. 
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